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ABSTRACT 

Using 7.5-year spectroscopic monitoring data of a sample of 17 Palomar-Green QSOs 
(PG QSOs) (z=0. 061-0. 371), we obtain the optical spectral slope for each object at all 
epochs by a power-law fit to the spectra in continuum bands. All of these 17 PG QSOs 
exhibit obvious spectral slope variability. Most of the 17 objects show anti-correlation 
between the spectral slope and the rest-frame 5100A continuum flux while five of 
them exist strong anti-correlation (correlation coefficient R larger than 0.5). For the 
ensemble of these 17 PG QSOs, a strong anti-correlation between the average spectral 
slope and the average rest-frame 5100A luminosity is found while no correlation is 
found between the spectral slope and the Eddington ratio. A median anti-correlation 
between spectral slope changes and continuum flux variations is also found which indi- 
cates a hardening of the spectrum during bright phases. Accretion disk (jet) instability 
models with other mechanisms associated with changes in the accretion processes are 
promising. 
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1 INTRODUCTION 

Variability phenomenon in active galactic nuclei (AGNs) 
provides a powerful constrain on their central engine. Var- 
ious classes of AGNs show variability on timescales from 
hours to decades and from X-ray to radio. The X-ray vari- 
ability in AGNs provides the strong constrain on the X- 
ray emission mechanics, the nucleus-size, the supermassive 
black hole, et al. (e.g. Mushotzky et al. 1993; Bian & Zhao 
2003). The large-amplitude and short-timescale optical vari- 
ability found in Blazars are likely due to the relativistic 
beam effects (e.g. Fan & Lin 2000; Vagnetti et al. 2003). For 
the small-amplitude and long-timescale optical variability in 
AGNs, its origin is still a question to debate and there are 
mainly four models: accretion disk instabilities (e.g. Rees 
1984; kawaguchi et al. 1998), supernovae explosions (e.g. 
Cid Fernandes et al. 1996; Aretxaga et al. 1997), gravita- 
tional lensing (e.g. Hawkins 1996), and star collisions (e.g. 
Torricelli-Ciamponi et al. 2000). It is still difficult to disen- 
tangle various models. 

For a long time the long-term multi-wavelength pho- 
tometric monitoring plays an important role in obtaining 
the optical variability of QSOs (see Table 1 in Giveon et al. 
1999). Many results on correlations between the variability 
and the luminosity, the redshift, the time lag, the rest-frame 
wavelength, spectral properties are obtained (e.g. Huang et 



al. 1990; Giallongo et al. 1991, Cimatti et al. 1993, Cristiani 
et al. 1996, Giveon et al. 1999, Trevese et al. 2002, Vanden 
Berk et al. 2004). Using the photometric data in various 
optical bands, the spectral variability can also be studied. 
Giveon et al.(1999) found that about half of the 42 Palomar- 
Green QSOs (PG QSOs) become bluer (harder) when they 
brighten (see also Trevese et al. 2001). Recently, using the 
photometric data from the Sloan Digital Sky Survey (SDSS), 
Vanden Berk et al.(2004) found quasars are systematically 
bluer when brighter at all redshifts for a sample of over 25000 
QSOs (see also de Vries et al. 2003). 

Although the photometric method have some advan- 
tages that the variability amplitude is accurate and that 
it takes less time comparing with spectroscopic method, it 
just monitors the flux variability in a few points in the wave- 
length which usually consists of many components includ- 
ing continuum and strong emission lines. So it is necessary 
to study the optical spectral variability using spectroscopic 
monitoring data. Recently, Wilhite et al.(2005) presented 
315 significantly varied quasars sample from multi-epoch 
spectroscopic observations of SDSS. They found that the av- 
erage difference spectrum (bright phase minus faint phase) is 
bluer than the average single-epoch quasar spectrum. How- 
ever, the epochs are just a few (~ 2) for these 315 SDSS 
QSOs. 

The AGNs spectral variability has focused on about 30 
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individual objects for a long-term of decades, with the pri- 
mary goals to obtain the structure of the broad line regions 
(BLRs) and the masses of their central supermassive black 
holes (Reverberation mapping methods, Peterson 1993). 
Kaspi et al. (2000) presented results from a spectrophoto- 
metrically monitoring program of a well-defined sample of 
28 PG QSOs. 17 of these 28 PG QSOs were observed with 
good time sampling (> 20 observing epochs). Good time 
sampling and long-term monitoring make these 17 QSOs the 
best objects in the studying of the spectral slope variability. 
In section 2, the sample, observations and data reduction are 
briefly described. The spectroscopic data analysis is given in 
section 3. Our discussions and results are presented in the 
last section. All of the cosmological calculations in this paper 
assume Ho — 70 km s _1 Mpc -1 , qo = 0.5. 



2 OBSERVATIONS AND DATA REDUCTION 

The 7.5-year spectra of 17 PG QSOs are available on the 
web 1 . The sample, observing technique, and reduction pro- 
cedure are briefly described here. The optical spectropho- 
tometric observations of 17 PG QSOs were done using 2.3 
m telescope at Steward Observatory (SO) and 1 m tele- 
scope at Wise Observatory (WO). The observations were 
performed between 1991 and 1998. The typical spectrum 
wavelength coverage at both observatories is from 4000 to 
8000A with a spectral resolution of about lOA. The redshifts 
of all these 17 QSOs are less than 0.4 and B < 16 mag. 
Total exposure times were usually 40 minutes at SO 2.3 m 
telescope and 2 hours at WO 1 m telescope. Spectrophotom- 
etry was obtained every 1-4 months for 7.5-year per object. 
Standard IRAF routines were used to perform spectroscopic 
data reduction. The consecutive quasar/star flux ratios were 
compared to test for systematic errors in the observations. 
Spectrophotometric calibration for each quasar was accom- 
plished excellently and the accuracies of order 1-2% can be 
achieved. Spectra were calibrated to an absolute flux scale 
using observations of spectrophotometric standard stars on 
one or more epochs. These 17 PG QSOs are primarily used 
to obtain BLRs size and the masses of their central super- 
massive black holes (Kaspi et al. 2000). The observed spec- 
tra have been corrected for Galactic extinction using Ab 
values from NED (See table 1 in Kaspi et al. 2000) 2 , as- 
suming an extinction curve with Rv — 3.1. Most of these 17 
QSOs are radio-quiet, and there are two radio-loud QSOs: 
PG 1226+023 (3C 273), PG 1704+608 with the radio loud- 
ness of 1621, 563, respectively (Nelson 2000). In Table 1, 
the object name, redshift, apparent B magnitude, number 
of spectrophotometric observing epochs are listed. 



3 DATA ANALYSIS 

Our goal is to investigate the long-term spectral variability 
of the PG sample. It is popularly accepted that we can use 

1 http: / /wise-obs. tau.ac.il/~shai/PG/ 

2 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 



the power law formulae, /„ oc v ~ a (/a oc A- 2+q ), to approx- 
imately fit the optical continuum spectrum of AGNs (e.g. 
Wilhite et al.2005). Some authors also used two power laws 
to model the continuum emission(e.g. Forster et al. 2001; 
Shang et al. 2005). The spectral slope usually changes with 
luminosity, which can naturally explain the relation between 
the variability and the wavelength. From the spectra of these 
17 PG QSOs during 7.5-year spectroscopic observations, we 
can obtain the spectral slope a (/„ oc v~ a ) by fitting a power 
law to the continuum spectrum, using spectral regions unaf- 
fected by other emission components. The usually used "con- 
tinuum windows" (at the rest-frame) known to be relatively 
free from strong emission lines are 3010-3040, 3240-3270, 
3790-3810, 4200-4230, 5080-5100, 5600-5630, 5970-6000, and 
6005-6035A(Forster et al. 2001, Vanden Berk et al. 2001). 
Here we directly used the continuum bands in the observer's 
frame defined by Kaspi et al. (2000) (See table 2 in Kaspi et 
al. 2000). The Balmer lines and other strong broad emission 
lines are excluded in these continuum bands. Most of these 
17 PG QSOs have 6 continuum bands. Two objects, PG 
1613+658 and PG 1700+518, have three continuum bands. 
From the best fit, we derive the spectral slope a, the con- 
tinuum flux at rest-frame 5100A and their errors for all ob- 
serving epochs per object. The average spectral slope and 
average flux density at rest-wavelength 5100A for each of 
these 17 PG QSOs are listed in Col. (5) and (6) in table 1. 
Here we mainly study the relation between a and continuum 
flux log / 5 ioo • 

In Fig. 1, we showed the spectral slope a as a func- 
tion of the rest-frame 5100A continuum flux log /5100 for 
each object. Considering the errors, it is obvious that all 
of these 17 PG QSOs exhibit spectral slope variability. We 
used the simple least-squares linear regression (Press et al. 
1992) to study the correlation per object. For each object, 
we obtained the correlation coefficients (R) and the prob- 
abilities (P) for rejecting the null hypothesis of no correla- 
tion. The fitting results for all spectral data per object are 
listed in Col. (7) and (8) in table 1. These best fittings are 
also plotted in each panel (Fig. 1, solid lines). From Fig. 
1 and table 1, we found that most of the objects showed 
negative correlation coefficients. Among the 15 QSOs which 
showed negative correlations, five objects PG 0026+129, PG 
0052+251, PG 0804+761, PG 1617+175, and PG 1704+608 
showed strong anti-correlations with R > 0.5. Two objects, 
PG 1229+204 and PG 1700+518, existed moderate oppo- 
site trends. One radio-loud QSO, PG 1704+608, showed the 
strongest anti-correlation with R — —0.88. While the other 
radio- loud QSO, PG 1226+023 (3C 273), showed weaker 
anti-correlation with R = —0.23. 

In order to give the ensemble relation for these 17 PG 
QSOs between the spectral slope and the luminosity, in the 
left panel in Fig. 2, we showed the average spectral slope 
a as a function of the average rest-frame 5100A luminos- 
ity log L5100 ■ We found that there existed a strong anti- 
correlation (R=-0.57, P=0.016). A simple least-squares lin- 
ear regression gives the relation, 

a = (-0.38 ± 0.07) log L5100 + (11-97 ± 2.08) (1) 

The best fitting is showed as solid line in Fig. 2. We 
also calculated the average Eddington ratio for each ob- 
ject, Lboi/ 'LEdd- The bolometric luminosity Ltoi is from lu- 
minosity in rest-frame 5100A: Lbol = 9ALa(5100A), and 
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L Edd = 1.26 x 10 38 M B h/M q ergs s~\ where the black hole 
mass is from Kaspi et al. (2000). In the right panel of Fig. 2, 
we plotted average a versus average Eddington ratio for each 
object. We found no correlation between them (R=0.064, 
P=0.081). 

Using the fitting slope of these 17 PG QSOs at all 
epochs, we correlated the spectral slope changes with the 
continuum flux variations. The changes of the spectral slope 
are calculated for each object between the slope of each ob- 
serving date and the average spectral slope. In Fig. 3, the 
overall 8a versus 51og/sioo is plotted. We found an anti- 
correlation between the spectral slope changes and the con- 
tinuum flux variations with R = —0.27 and P < 10 -4 . A 
simple least-squares linear regression gives, 

8a = (-0.61 ± 0.01) x 8 log /5100 - (5.91 ± 0.52) x 10~ 3 (2) 

The best fitting is showed as the solid line in Fig. 3. 



4 DISCUSSIONS AND RESULTS 

These 17 PG QSOs data span 7.5 years with more than 20 
epochs and are used in the reverberation mapping method to 
determine their central supermassive black hole masses (e.g. 
Peterson 1993; Kaspi et la. 2000). The time-span is enough 
to determine the variability of the optical spectral slope if 
any. For these high-luminosity PG QSOs, the contribution 
from the host galaxies is negligible. Here we used a power 
law function to fit data in the continuum bands suggested by 
kaspi et al. (2000). And then we obtained the spectral slope 
and the flux /5100 in the rest-frame. Their errors are from 
our fitting. The results would possibly depend on the selec- 
tion of the continuum bands. In order to be consistent with 
the continuum analysis of Kapsi et al. (2000), We adopted 
continuum bands suggested by kaspi et al. (2000) instead of 
usual continuum windows(Forster et al. 2001, Vanden Berk 
et al. 2001). In the future we will use the multi-component 
fitting method to avoid the selection effect of the continuum 
windows (Bian, Yuan & Zhao 2005). 

4.1 Correlation for individual objects 

Giallongo et al.(1991) investigated the long-term optical 
variability of quasars from photometric observations and 
found a positive correlation between variability and redshift. 
That QSOs at higher redshifts show larger variability since 
they are observed at a higher rest-frame frequency where the 
variability is stronger was indicated. This suggests a hard- 
ening in the bright phase of quasars. Statistical evidence for 
this trend were presented by follow-up studies (e.g. Cristiani 
et al. 1996, Di Clemente et al. 1996). Giveon et al.(1999) pre- 
sented results from 7 yr photometric monitoring program of 
a well-defined, optically selected sample of 42 PG QSOs. 
The spectra of about half of the QSOs in their sample be- 
came harder when they brighten. The 42 objects in the op- 
tical sample included all 17 objects discussed in the present 
paper. Here on the basis of 7.5-year spectroscopic monitor- 
ing data of these 17 PG QSOs, we correlated the spectral 
slope with the rest-frame 5100^4 continuum flux. The re- 
sults showed that, among these 17 PG QSOs, five objects 
displayed strong anti-correlation with R > 0.5 and there are 
nine objects with R > 0.3 (See Fig. 1 and table 1). The 



anti-correlation implied that QSOs become bluer as they 
brighten. At the same time, other objects showed weaker 
correlation between a and /5100 (See Fig. 1 and table 1). For 
PG 1700+518 and PG 1229+204, the correlation is possibly 
positive. The Fe II multiples are strong in PG 1700+518 
(see Fig. 1 in Kaspi 2000). The strong Fe II emission may 
contaminate the selected continuum bands and the studying 
of the spectral slope variability for this object is probably 
influenced. The spectra of these 17 PG QSOs are from 2.3 
m telescope at Steward Observatory and 1 m telescope at 
Wise Observatory. For each object, a large fraction of the 
spectral data are from the Wise Observatory. Considering 
only the spectra from the WO, a strong anti-correlation be- 
tween the spectral slope and the continuum flux is found 
in PG 1229+204 (R = -0.80, P < 10~ 4 ) (see Fig. 1). The 
least-squares linear fitting using the spectral data only from 
the WO are also showed for all the other PG QSOs in our 
sample (Fig. 1, dotted lines). The result for this fitting are 
listed in Col. (9) and (10) in table 1. Using the spectra only 
from the WO, all but the object PG 1700+518 showed anti- 
correlation between the spectral slope and the continuum 
flux. Stronger anti-correlation are obtained for some of the 
objects (e.g. PG 1426+015, PG 1613+658, PG 2130+099). 

4.2 Correlation for the ensemble of 17 objects 

In order to clarify the slope variability, some authors dis- 
cussed the relation between the spectral energy distribution 
(SED) changes and the continuum flux variations. Using the 
variation between two epochs in the photographic U, Bj, F, 
and N bands, Trevese et al.(2001) correlated the spectral 
slope changes Aa with brightness variations A log /./ for a 
complete magnitude-limited sample of faint quasars in SA 
57 and detected an average increase of the spectral slope 
(note that their a is defined by /„ tx v a ) for increasing con- 
tinuum flux, indicating hardening of the spectrum in the 
bright phases. Trevese et al.(2002) performed an analysis 
of B and R observing data of 42 PG quasars from Giveon 
et al.(1999). They showed in their results that the average 
spectral slope of each QSO tended to be larger for brighter 
objects (the a is defined by /„ cx v a ). Using 7.5-year spec- 
troscopic observing data, we discussed the relation between 
Aa and Alog/5100 for the ensemble of these 17 PG QSOs, 
a median anti-correlation (R = —0.27, P < 10 -4 ) between 
them is found (see Fig. 3), indicating a hardening of the 
spectrum during bright phases which is consistent with the 
result of Trevese et al.(2001). We also found a strong anti- 
correlation between the average spectral slope and the aver- 
age rest- frame 5100A luminosity (See the left panel in Fig. 2) 
which is qualitatively consistent with the inter-QSO result 
of Trevese et al. (2002). 

4.3 Model behind optical spectral variability 

The physical mechanisms behind the optical variability of 
AGNs is largely unknown. Some models have been proposed 
to explain the AGNs optical variability properties. Current 
models can be classified mainly into three groups: accretion 
disk instabilities, discrete-event or Poissonian processes, and 
gravitational microlensing (e.g. Hawkins 1996, Aretxaga et 
al. 1997, kawaguchi et al. 1998, Kong et al. 2004). The opti- 
cal spectral variability can provide clues to disentangle these 
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models. Using the optical spectra of NGC 5548 between 1989 
and 2001, Kong et al.(2004) found a strong correlation be- 
tween a and the luminosity at rest-frame 5100A. They used 
the global variance of the accretion rates or/and the vari- 
ance of the inner radius of the accretion disk to account for 
the strong relation. It can explain well the anti-correlation 
between the optical spectral index and the continuum lumi- 
nosity observed in NGC 5548. Here in our PG QSOs sam- 
ple, this disk instability model can also explain the anti- 
correlation which most of the objects showed between the 
spectral slope a and the rest-frame 5100A continuum flux for 
individual QSOs. However, for the ensemble of our PG QSOs 
sample, the global variance of the accretion rates would not 
necessarily account for the global slope variance (See the 
right panel in Fig. 2). Vanden Beck et al. (2004) studied 
the ensemble photometric variability for about 25000 SDSS 
QSOs, that quasars are systematically bluer when brighter 
at all redshifts were found. Their result also seemed to dis- 
favor gravitational microlensing and generic Poissonian pro- 
cesses as the primary source of quasar variability. QSOs are 
widely believed to be powered by the accretion disk onto 
a supermassive black hole (e.g. Rees 1984). It is natural to 
consider that the QSOs variability is due to some mecha- 
nisms associated with changes in the accretion processes. 
Kawaguchi et al. (1998) presented a very simple cellular- 
automaton model for disk instability. They gave the slope 
of structure function between 0.41 to 0.49, which is incon- 
sistent with SDSS results (Vanden Beck et al. 2004). It is 
likely due to the complexity of possible accretion disk (or 
jet) instability models, which prevented more quantitative 
predictions. 



5 CONCLUSION 

Using 7.5-year spectroscopic monitoring data of a sample of 
17 PG QSOs, we study the optical spectral slope variance. 
The main conclusions can be summarized as follows: 

• Using the continuum bands suggested by Kaspi et al. 
(2000), we found that, in 7.5-years long-term observation, 
all 17 PG QSOs showed obvious optical spectral slope vari- 
ability. 

• Most of these 17 PG QSOs showed anti-correlation be- 
tween the spectral slope a and the rest-frame 5100A con- 
tinuum flux log /5100 while five of them showed strong anti- 
correlation between a and log/5100 (R > 0.5). 

• For the ensemble of these 17 PG QSOs, a strong anti- 
correlation (R — —0.57) between the average spectral slope 
a and the average rest-frame 5100A luminosity logZ/5100 is 
found while a median anti-correlation (R = —0.27) is found 
between spectral slope changes 5a and continuum flux varia- 
tions 5 log /5100 indicating a hardening of the spectrum dur- 
ing bright phases. 

• The disk instability model can qualitatively explain the 
anti-correlation between the spectral slope a and the contin- 
uum flux log /5100 for individual QSOs. However, the global 
variance of the accretion rates would not necessarily account 
for the global slope variance for the ensemble of our PG 
QSOs sample. Accretion disk (jet) instability models with 
other mechanisms associated with changes in the accretion 
processes are promising. 
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Table 1. Sample properties. Col. l:The name, Col. 2: rcdshift, Col. 3: apparent B magnitude, Col. 4: number of spcctrophotomctric 
observing epochs, Col. 5: the average optical spectral slope, Col. 6: average continuum flux density at rest wavelength 5100A in units of 
10 — 16 ergs s~ 1 cm~ 2 A -1 , Col. 7-8: correlation coefficient and probability for all spectra per object, Col. 9-10: correlation coefficient and 
probability for the spectra only from the Wise Observatory. 
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Figure 1. The spectral slope a versus the rest-frame 510CL4 continuum flux log/5100 (in units of 10~ 16 ergs s _1 cm _2 A~ 1 ) for each of 
these 17 PG QSOs. The spectral data from the Steward Observatory are solid rectangles, and the spectral data from the Wise Observatory 
are open rectangles, respectively. Solid lines are the best fit to all spectral data for each object while dotted lines are the best fit to the 
spectral data only from the Wise Observatory. 




Figure 1. Continued. 
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Figure 2. Left: the average spectral slope a versus the average rest-frame 5100A luminosity logZ/5100 in units of ergs s~ 1 Hz~ 1 for 17 
PG QSOs. Right: the average spectral slope a versus the average Eddington ratio for 17 PG QSOs. The solid lines are our best fits. The 
correlation coefficient R and the probability for rejecting the null hypothesis of no correlation are showed in the right figure corners. 




Figure 3. The spectral slope changes 6a versus the rest-frame 510CL4 continuum flux variations (Jlog/sioo- The solid line is our best fit. 



